Magnetic field induced inequivalent vortex zero modes in strained graphene 
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Zero energy states in the Dirac spectrum with U(l) symmetric massive vortices of various under- 
lying insulating orders in strained graphene are constructed in the presence of the magnetic field. 
An easy plane vortex of antiferromagnet and quantum spin Hall orders host two zero energy states, 
however, with two different length scales. Such inequivalent zero modes can lead to oscillatory 
charge and magnetization, and their usual quantizations get restored only far from the vortex core. 
Otherwise, these zero modes can be delocalized from each other by tuning the mutual strength of 
two fields. One can, therefore, effectively bind a single zero mode in the vortex core. A possible 
experimental set up to capture signature of this theory in real graphene as well as in optical hon- 
eycomb lattices is mentioned. Generalization of this scenario with underlying topological defects of 
Kekule superconductors can localize a single Majorana mode in the vicinity of the defect-core. 

PACS numbers: 71.10.Pm, 71. 10. Li, 05.30.Fk, 74.20.Rp 
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I. INTRODUCTION 



The existence of states at precise zero energy in the 
spectrum of a relativistically invariant Dirac equation 
with underlying topological defects [i| , recently attracted 
ample interest, following the successful fabrication of 
graphene. The extent of the spectrum over the positive 
and the negative energies provides robustness to such 
states at the middle of the spectrum against weak 
local perturbations. The number of zero energy states, 
however, depends only on the universal properties: total 
magnetic flux enclosed by the system 0], kink of the 
mass ^3;], and the vorticity of a mass vortex For 
example, a half-filled Landau level at zero energy, leads 
to the formation of Hall plateaus at fillings v = ±2 
The existence of zero energy modes can lead to 
fractionalization of quantum numbers, e.g., charge 
Majorana modes in the core of the superconducting 
vortex or half-vortex and additional competing 
order parameters in the core of the vortex Q. 
The Kosterlitz-Thouless (KT) scaling of longitudinal 
resistivity {Rxx) in neutral graphene strongly suggests 
the possibility of vortex excitations in graphene subject 
to the magnetic fields [l0|. Vortex-zero modes in the 
Dirac spectrum can exist in the presence of either 
real or pseudo magnetic fields The nature of the 

electronic ground state in the presence of fictitious 
gauge field lately gained attention following the exper- 
imental realization of the strain induced fields [12h14| . 
In the presence of uniform real {B) and pseudo (6) 
magnetic fields, quantized Hall conductivity in graphene 
is expected to discern plateaus at all integer values of 
e^//i, however, at noninteger fillings, when the former 
is kept stronger. The plateau of Hall conductivity 
with Uxy — e^/ft., can be formed by developing a 
ferrimagnet order, due to distinct degeneracies of two 
inter-penetrating relativistic Landau levels living in 
the vicinity of two inequivalent Dirac points. Other- 
wise, their degeneracies are proportional to the effective 



gauge fields near the Dirac points, B±h, respectively tl5| . 

A fundamental question is, therefore, raised: Does 
the vortex like defect have any zero energy mode when 
graphene is subject to both the real and the pseudo 
magnetic fields ? Deposition of graphene on a metallic 
substrate at relatively high temperatures, followed 
by cooling, produces a strain induced pseudo (axial) 
magnetic field due to the mismatch in compressibility 
with the substrate [12|. This system can then be 
placed in a real magnetic field. Recently, a hexagonal 
honeycomb lattice has been realized in an optical 
lattice system (0LS)[T6'|. One can, therefore, study 
the proposed scenario in honeycomb OLS as well. In 
the OLS, a magnetic field is introduced via synthetic 
gauge potentials, and its strength can be varied over a 
wide range On the other hand, a pseudo magnetic 
field can arise from specific modulations of the nearest- 
neighbor hopping (NNH) amplitudes only [l3l. [l l [Tsl fioj . 

In graphene possible candidates with requisite U(l) 
symmetry are antiferromagnet (AF), Qquantum spin 
Hall orders projected onto the easy plane and the Kekule 
bond density wave (KBDW) . The former two can acquire 
an easy plane in the magnetic field due to the Zeeman 
coupling [l^l. Otherwise, AF and QSH phases are 
energetically favored by strong onsite [22| and second- 
neighbor repulsions, respectively [IJ. KBDW can be 
stabilized by electron-phonon interactions (23 |. However, 
a QSH phase can be realized even at sufficiently weak 
next-neighbor repulsion when a fictitious magnetic field 
penetrates graphene [ij, [l^. Immaterial of the vortex 
nature of Dirac quasi-particles we find that there are 
always two zero energy modes. However, categorically 
two different situations arise: When the underlying order 
parameters do not couple two inequivalent Dirac points, 
e.g., AF and QSH orders, zero energy modes have two 
different characteristic lengths. With underlying KBDW 
order, two zero energy modes differ from the ones in the 
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absence of pseudo fields only by an overall factor, hence, 
enjoy a unique length scale. In the former situation, we 
show that the zero modes can be delocalized from each 
other arbitrarily by tuning the relative strength of the 
fields only. The present mechanism can, therefore, allow 
one to localize a single zero mode near, while pushing 
the other one far from, the vortex-core. The charge and 
the magnetization then exhibit oscillatory behavior, 
while restoring their usual quantization only at large 
distances. A previous study shows that in the absence of 
gauge fields the zero energy modes have the same decay 
leng ths, hence, exact quantization is restored everywhere 
[25| . Additionally, if there exists a half-skyrmion of 
AF order, however weak, which otherwise splits the 
zero modes if uniform, maximally lowers the energy 
by filling the state with shorter length scale. The zero 
energy states can be computed exactly if the fields are 
uniform, although they always exist for any arbitrary 
flux profile and/or shape of the vortex mass. Tuning the 
strength of the synthetic gauge fields and/or modulated 
NNH, thereby of pseudo magnetic fields, one can grasp 
peculiarities of these zero modes on honeycomb OLSs. 

The rest of paper is organized as follows. In the 
next section, we write down the free Dirac Hamiltonian 
subject to both real and pseudo magnetic fields and 
mention the single-particle Landau-level spectrum. 
Existence of zero energy modes in the presence of an 
underlying vortex with easy-plane AF order is discussed 
in Sec. III. Robustness of the zero modes against any 
modulation of the field and/or the vortex mass profile is 
presented in Sec. IV. Section V is devoted to study the 
internal structure of the two-dimensional zero energy 
manifold and various competing orders in the core of 
the vortex. Mid-gap states with underlying vortices of 
easy-plane QSH order and KBDW order are respectively 
considered in Sees. VI and VII. Concluding remarks and 
a discussion on related issues are presented in Sec. VIII. 



the relativistically invariant Hamiltonian takes the form 
Hu = I2 ® ijoJiPi- The four-component mutually 
anticommuting Hermitian 7-matrices belong to the 
"graphene representation", 70 = ^2 (H) , 71 = CT3 (g) cr2 , 
72 = /2<8)o-i, 73 = cri(Ki(T2, and 75 = a2®cT2- Here, (/2,<?) 
are the standard Pauli matrices. The chiral ?7c(4) sym- 
metry of Hd is generated by {h,$) ® ih,"/3,l5in3l5) 
[2^ . Three generators ofrotation of electron spin are 
S ^ I4. 

The free Dirac Hamiltonian in the presence of both 
real (Ai) and pseudo (a^) gauge potentials reads [11] 

Hd [A, a] = /2 Z7o7i [Pi - M - ilsl^ai) ■ (2) 

The real (pseudo) magnetic field breaks (preserves) 
the time reversal symmetry (TRS), represented by 
It ~ (J2K ® ijijc^K in "graphene representation" and 
K is the complex conjugate [13, [13] j but preserves 
(breaks) chiral symmetry (CS). Separately, both real 
and pseudo magnetic fields quench the linear spectrum 
of Dirac quasi-particles into a set of Landau levels at 
well separated energies ±y^2nB{b), with n ~ 0,1,2,---, 
and degeneracies J727ri?(&), where ft is the area of the 
sample. However, the states in zeroth Landau level 
near two Dirac points reside on the complimentary and 
same sub-lattices, respectively. The distinct natures of 
the zeroth Landau level can lead to different broken 
symmetry phases near and at the charge-neutrality 
point [13, [ij]. On the other hand, in the presence of 
a real as well as a pseudo magnetic field, the spectrum 
of the Dirac quasi-particle is composed of two inequiv- 
ale nt interpen etrating sets of Landau levels at energies 
±y/2n]B±b\, with degeneracies Qn\B ± b\, respectively, 
residing in the vicinity of two Dirac points '29']. Next, 
we consider vortex defects of various underlying U(l) 
symmetric order-parameters in graphene, subject to real 
and pseudo magnetic fields and study the zero energy 
modes bound to them. 



II. FREE HAMILTONIAN WITH REAL AND 
PSEUDO MAGNETIC FIELD 

To describe the massless, chiral Dirac fermionic 
excitations around the two inequivalent corners of 
the Brillouin zone, suitably chosen at Ki — 
-K2 = (l,l/\/3)(27r/a\/3) ^26^, we construct an 
eight-component spinor '^{x) 



{^^{x),-^^{x))^ , with 



{u\^{x) , (f ) , {x) , vl^ [x) ) 



stands for electrons spin projection along the z— direction 
and 



YUx) 



dp 



(2™) 



'p-^YUi-^r'K+p). (1) 



Y = u,v are the electron annihilation operators on 
sublattices A and B, respectively. In this representation. 



III. VORTEX WITH EASY-PLANE 
ANTI-FERROMAGNET ORDER 

The onsite Hubbard repulsion is the strongest interac- 
tion in graphene [13] and favors AF order at neutral filling 
if sufficiently strong (32| . Hence, we consider a uniform 
background of electron density and staggered magnetiza- 
tion. Therefore, after keeping the only relevant term, the 
onsite Hubbard interaction reads as 



Hu = -^[f{A)-f{A + b)]\ 



(3) 



Here, f{X) = cl{X)(Jaa'Ca'{X) is the average magne- 
tization on sublattice X = A,B, defined in terms of 
= Ua{A),Va{B), respectively. The effective single- 
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particle Hamiltonian in a fixed Neel background is 

HN[A,a]=HD[A,a]-[N{x)-a)®-to, (4) 

where N = {f{A)~f{A + b)) 7^ [H IH . In the absence 
of gauge fields, a constant Neel order, N{x) = N, leads 
to a gapped spectrum for the relativistic quasiparticles, 
E = ± -\/fc^"^r~pVp, and the order parameter (TV) breaks 
the CS [2^. In the presence of real magnetic fields, a 
finite Neel order shifts the Landau levels at finite energies 
to ±y/2nB + N"^ with degeneracies 27ri?, whereas, the 
zeroth Landau level gets splits to ±.\N\ with degeneracy 
ttS, per unit area [21|. Next, we assume N^^^x) ~ 
and Ni{x) + iN2[x) = |A(r)|e^P'^, p = ±1: vortex 
configuration [s^ ]. Here, we use a rotationally invariant 
representation of spinor '^cr{x). Thus the choice of the 
easy plane is arbitrary and does not affect the outcomes. 

In the absence of gauge fields, there exist two zero en- 
ergy modes in the spectrum of i?Ar[0, 0] d^]. Neither 
H£)[A,a\ nor 70 mixes two inequivalent Dirac points, 
hence, H^lA^a] is block diagonal in the valley index. 
One can, therefore, cast H^lA^a] as Hi ® H2, by ex- 
changing the second and third 2x2 blocks. Both Hi and 
H2 however, are unitarily equivalent to a generic Dirac 
Hamiltonian with mass-vortex in the presence of different 
effective magnetic fields. 



H = HoiA') + Ni{x)i-fol3 + iV2(f)i7o75, 



(5) 



-id, - Al). 
ia2 and A*j 



Specifically, Hi ~ 
= Ai -\- aj , whereas 



where Hq{A^) = i7o7i( 
uIhUi, with Ui= 
H2 = UIHU2, with U2 = ia2 e I2 and Aj = A, ~ a, 0. 
Therefore, the effective single-particle Hamiltonian with 
a vortex for the easy plane of the Neel vector is equiva- 
lent to two copies of the Dirac Hamiltonian with a twisted 
mass (vortex), introduced by Jackiw and Rossi [Tl| and re- 
cently studied in the context of Kekule orders in g raphene 
0, but subject to effective magnetic fields [IJ. These 
two copies, however, experience different effective fields. 
For example, choosing XA{r) — Br/2 and Xa(^) = ^ 
so that A{r) = XA{r){y,~x) and a{r) = Xa{r){y,-x) 
yields uniform magnetic fields B ± b for -^1(2)- Since 
the pseudo magnetic fields can be produced by deliber- 
ately buckling the graphene flake, it is likely to be non- 
uniform, in general. Nevertheless, a previous study al- 
ways guarantees the existence of a zero energy modes in 
the spectrum of Hi and H2, irrespective of the modula- 
tion of the fields and/or vortex profile pjj. Therefore, 
Hn[A, a] = Hi Q) H2 always hosts two zero modes. En- 
couraged by a recent experiment [l2j, we can take the 
axial field to be uniform. The characteristic lengths of 
two zero modes depend on the effective uniform magnetic 
fields. Therefore the zero energy sub-space is composed 
of two modes with different length scales. Let us consider 
an extreme limit: B zL b ^ ^('') = ^0 (constant). To 
the leading order, one can write the localization lengths 



of the zero modes as 



o(. 



A2 
^0 



VB±b ^B±b 



(6) 



This result physically makes sense, since with a weak 
order parameter, to the leading order, magnetic fields 
determine the characteristic length scales. Next, we 
consider a prototype vortex profile, |A(r)| — rAp/i? 
n (const 
[ll| an( 

modes explicitly and thereby L±. 



when r < R and Aq (constant) when r > R and the 
fields to be uniform 11 il l and compute the zero energy 



IV. ZERO MODES IN UNIFORM FIELDS 

Before we arrive at the zero modes for the easy-plane 
Neel vortex, it is worth studying the zero modes with the 
Kekule mass vortex in the presence of an effective real 
magnetic field. Even though this problem has been stud- 
ied previously [ll|, it is worth reviewing the solution for 
the sake of completeness of the discussion. The equations 
for the zero energy mode, '^^oix) = {u+,v+,U-,V-)^ , of 
the massive Dirac Hamiltonian, 

H = i7o7i {Pi - Al) -\- I A(r)| (17073 cos -1- 17075 sin 0) 

(7) 

with a mass vortex and effective magnetic field reads as 



drqir) 
drp{r) = 



-\A{r)\p{r) ^xHr)q{r), 
-\A{r)\q{r)+x\r)p{r). 



(8) 
(9) 



Here we choose a symmetric gauge A^^r) — %*(''') (^j 
and redefine the spinor components as u+(r) = y/i p(r) 
and u_(r) = a/— « q{r), with v± — 0. The other possi- 
bility with u± = 0, does not lead to a normalizable zero 
energy state as in the absence of fields [1, HI] . For exam- 
ple, assuming a finite field at the origin, i.e., xi'"') oc r, 
yields v± cx 1/r near the origin. Therefore, v± cannot 
be normalized. To find the zero modes, let us define a 
time like variable, dt = |A(r)| dr, as in Ref. [ll|. One 
can then identify the above two equations as the Hamil- 
ton's equations, previously argued in Ref. [ll|, for the 
canonical momentum (p(r)) and position (g(r)). The cor- 
responding time varying classical Lagrangian, 



1 d 



2 '^-+2^[^*W'^^]' 



describes the motion of a classical particle in a repulsive 
harmonic potential (upto a total time derivative), 
C{t) = 1 + f\ty - fit), where f\t) = x\t)/\A{t)\. 
In the presence of magnetic fields the potential gets 
steeper with time, otherwise, it remains static. In the 
mechanical analogy, one may understand the zero energy 
mode as the trajectory of a particle starting far from the 
origin with the right amount of energy, so that it reaches 
the origin with zero velocity at infinite time. Since the 
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particle with sufficiently large (small) initial energy flies 
(goes back) to negative (positive) infinity, continuity of 
solutions, therefore guarantees the existence of the zero 
energy mode in the spectrum H in Eq. [7]. 

The vortex Hamiltonian with the easy plane compo- 
nents of the Neel order, subject to magnetic fields, is 
shown to be unitarily equivalent to two copies of H in Eq. 
[7], namely, Hi ® H2. However, x^i^) = XA{r) + Xa{r) 
for Hi and x*(^) = XA{r) - Xa{r) for H2. Here, x^(r) 
and Xaif) are the vector potentials of the real and 
the pseudo magnetic fields, respectively. The above 
mentioned mechanical analogy of the zero modes ensures 
the existence of two zero ener gy states, irrespective of 



the field and the vortex profile [11 1 



Next, we consider %*('') ~ {B ± h)r/2 for i?i(2)j re- 
spectively. Therefore, the effective fields read as i?e// — 
{dr + l/r) x*(r), yielding Beff ^ B + b = B+ ioi Hi and 
B — b = B- for H2 ■ Let us take the vortex configuration 
as in Ref. [ll|. |A(r)| = Aqv/R for r < i? (vortex-core) 
and Aq for r > R. For r < R, the components of two 
zero modes, say ^'i^o and ^'2,07 are 



;(r) = C±e'"/^i + C±e-'-^/^i, (11) 
and canonical momentum p{r) = dC{q,q,t) / dq is 
-2Rr 



plir) 



AoLl 



(12) 



Here, ± corresponds to the zero eigenstates of ^?^i(2)j re- 
spectively. Cf^ are constants, and L± are the character- 



istic lengths of the zero modes. 



1 fAi , iB±br 



r-4 _ - — U 



Whereas, far from the core of the vortex (r > R) 
(B±b).^ f3\S\ + l~S ^5_{B±b)i 



(13) 



g^(r) ^Cse —p 



A\S\ '2' 



(14) 



when S = B/2AI < 1. On the other hand, for 5 > 1 
>, , ^ _(B±^ f2S+l 3 {B±by\ . 



AS ' 2' 



where F is the hypergeometric function [ll|, |3l| 
r > R, dt = Agdr, yielding 



(1/Ao) {dq^{r)/dr). 



For 



(16) 



Continuity implies that the value and the first derivative 
of the solutions in two regions must match at r = i?, 
which eliminate two out of the three constants, whereas, 
the remaining one is fixed by the normalization condi- 
tion. Equation (IT5|) reflects that two zero energy modes 



have different characteristic lengths. Otherwise two zero 
modes of Hn[A, a] are 

vPi^o = C//(p+,0,(z+,0)^(S+), (17) 
*2.o = C/l(P-,0,g-,0)^(B_), (18) 

where Ui = 12® icr2 and U2 = ia2 ® h- Note VPi^o 
and ^'2,0 are the functions of the effective fields B + b 
and B — b, respectively. In the following discussion, 
we show how the existence of two length scales yields 
various peculiarities. However, the following discussion 
is insensitive to the exact form oi L±. 



V. INTERNAL STRUCTURE OF THE ZERO 
MODES 

In the presence of real and pseudo magnetic fields, the 
zero energy eigenstates of Hi and H2 are 

*i,o ("It, vit, uii, ui;)^ = Vi (p+(r), 0, 0, -iq+(r))^ , 

(19) 

and 

*2.o = {u2^,V2^,U2i,V2i)^ = (0, p_ (r) , (r) , 0)^ , 

(20) 

respectively. q±ir) and p±{r) have been computed 
exactly in uniform fields with the above mentioned 
profile of the vortex mass. Typically, the characteristic 
length for ^i^ is smaller than that for ^2,0i for example 
L± in Eq. (|13l) . Consequently, the former is squeezed 
into the vortex-core. 

The magnetization within the zero energy sub-space is 



Hz = X (era ® ^2) © (crs^h), 



(21) 



proportional to the Zeeman term. X = g B{x) and 
g fa 2 for electrons in graphene. In the absence of 
fields, p±{r) = q±{r) = exp[— |A(z)| dz]. Each of 
the states then possesses zero magnetization everywhere 
in the space [2^. In the presence of fields, the radial 
dependences of q{r) and p(r) are different, in general. 
Therefore, the zero modes can exhibit finite but oscilla- 
tory magnetization, shown in the inset of Fig. [T] The 
overall net magnetization is still zero for each state and 
zero energy sub-space remains unperturbed. Such a lo- 
cal magnetic moment can be probed by magnetic force 
microscope (MFM) measurements. Replacing the vor- 
tex by an anti- vortex simply exchanges the role of u and 
V components. For b > B, ^'i,o remains invariant, but 
the role of q-{r) and p-{r) is exchanged in ^2,0- This 
limit seems quite conceivable, since, in the experiment, 
b ^ 350 T [l2| , but the highest laboratory magnetic field 
is - 45 T. 

Two zero energy modes ^i = (^'i.o,0) and 
^2 = (0, ^'2.0) comprise a two dimensional basis 
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L, 



FIG. 1: Energy {E/N3) as a function of magnetic length 
{Lb /a), witli B{r) = B and N^ir) = N3 for r < R, zero 
elsewhere. The vortex configuration is A(r) — Ao for r > R, 
otherwise zero. We set R = 5 A and Ao — 30 eV. Inset: Ab- 
solute magnetization QM/msl) as a function distance (in A) 
from the center, for Lb ~ 40A ~ 

in the zero energy manifold (Ho). Any operator 
that commutes or anticommutes with Hp; [A, a] leaves 
that space invariant. When a = 0, there are 
four matrices falling in the second category, namely 
(0-3 <8) 70 , -^2 <8i 47073 , -^2 fX" «7o75 , 0-3 'S> i7i 72 ) , which to- 
gether, close a Cl{3) x U{1) algebra [1,|33|. 0-3(8)70 is the 
z-component of the Neel order and the remaining mem- 
bers of the Cl{3) group correspond to the KBDW. The 
U(l) part is generated by (T3 (g) 17172, representing the 
third component of the spin-triplet version of the time- 
reversal symmetry- breaking order . When the pseudo 
magnetic field is finite, the KBDW orders no longer anti- 
commutc with HD[A,a]. Nevertheless, the expectations 
value of all four orders come only from "Ho- The expec- 
tation value of any physical observable is 

= M E - E ) ^^Um^Eix), (22) 

\occup empty/ 

where M is the traceless matrix and {^£;(x)} is the 
set of eigenstates, with eigenvalues i? of a generic 
Hamiltonian. The existence of a unitary matrix T, 
which commutes with M, while anticommuting with 
the Hamiltonian, restricts the above sum within Hq 
[25I [36j . Choosing T = (73 70, which anticommutes 
with H]\[[A,a], one finds that the expectation value of 
either M = 0-3 (8) 70 or 0-3 (g) 27172 is solely determined 
by ^i,o's. For M = I2 x 27075 or I2 x 17075 or any linear 
combination, one can choose T — <^ «7i72- Therefore, 
the AF and QSH as well as the KBDW orders acquire 
their expectation values only from the zero energy 
sub-space. 

When the Neel vector is tilted out of the easy plane 
{N3 ^ 0), the zero energy manifold gets split, since 

A^3 (*Io (f^3 «) 70) *j,o) - ±N3S,j [25]. Therefore, at 
half-filling, a finite N3 keeps only one of the zero energy 
states occupied while leaving the other one empty. On 
the other hand, if the Neel order forms a half-skyrmion, 
^i^o's are no longer degenerate. If |A^3(r*)| <^ Ao, 
one can neglect the mixing of zero modes with the 
rest of the spectrum. The energy is then maximally 
lowered by filling the state with smaller characteristic 



length (L+). The physical reason: N^i^x) is finite in 
the core of the vortex and smoothly vanishes toward 
the boundary. Out of the two states, the one with the 
smaller length scale experiences a larger overlap with the 
out-of-plane component of the Neel order and, therefore, 
being filled lowers the energy maximally. The effective 
magnetic fields for '^1,0 and '^2.0 are B + b and \B — b\, 
respectively. The characteristic length is smaller for 
^'i^o- Hence, state ^'i.o, is occupied in the presence of a 
weak half-skyrmion or meron configuration of N. This 
feature is independent of specific configuration of field 
or vortex or N^. Nevertheless, energy as a function of 
magnetic length, for a chosen vortex profile, N3 and the 
fields is shown in Fig. [1] 



VI. EASY-PLANE QUANTUM SPIN HALL 
VORTEX 

If, on the other hand, repulsion among the fermions 
living on the second-neighbor sites of the honeycomb lat- 
tice is sufficiently strong, a QSH order {C — {J^a '^ti'^^ 
47172] ^(t) 7^ 0) can be realized in graphene. It breaks the 
TRS for each spin component and corresponds to circu- 
lating currents among the sites of the same sub-lattice 
[23|. However, the QSH phase can be realized even at in- 
finitesimal next-nearest-neighbor interactions in the pres- 
ence of finite pseudo magnetic fields jisl-fisl [Tsj . The 
single-particle Hamiltonian with an underlying QSH or- 
der reads as 

He [A, a] = Hd [A, a] - (c (f) • a) ® 271 72. (23) 

Alike AF, the QSH order is also block diagonal in the val- 
ley index and can be projected onto the easy plane by a 
finite Zeeman coupling [3J]. Otherwise, Ni — > {—ly^^d 
for Hj with i,j = 1,2 [1^. The two zero energy 
modes are the following: ^'i.o remains unchanged, while 
p_(r) o q-{r) in ^'2,0- Within Hq, the members of the 
Cl{3) algebra are (0-3 70, crs <g ^7073, 0-3 g) 27075 )• The 
last two entries corresponds to the z-component of the 
spin-triplet KBDW. The U(l) component is formed by 
C3 €5 27172- A finite third component of the QSH order 
(C3) can be developed in the core by making the vortex 
charged, since it acts like an identity operator in Hq- 
An out-of-plane component of the Neel order splits Hq. 
Therefore, onsite Hubbard repulsion (U) can lift the 
degeneracy of the zero energy sub-space by developing a 
finite A'3 in the core of the vortex. 

It is worth confirming that all four mass-orders, 
constituting the Cl{3) x U{1) algebra, acquire their 
expectation values from Hq, even though some of the 
members do not anti-commute with the Hamiltonian 
upon imposing an axial magnetic field. If we select 
M = (73(8)70 or (73 (g) 27172, one can choose T = 0-3 (g) 70. 
And for AI — <Si 27073 or (73 (g) 27075, one can again 
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choose T = 0-3 (g) 17172, same as before. 



VII. KEKULE VORTEX 

In contrast to the AF and QSH orders, the KBDW cou- 
ples two inequivalent Dirac points and breaks the trans- 
lational symmetry of the lattice into the Kekule pattern 
[1,0|- the presence of gauge fields, the effective single- 
particle Hamiltonian Hk[A, a] with an underlying vortex 
of KBDW order is 

HD[A,a]-mrl2(^ijo{73Cos(j> + j5Sm(l)) = UHk[A,0]U, 

(24) 

where U = exp[— p(a;)(/2 (8) 7o)]- nT-r counts the radial 
variation of the vortex profile. I2 fXi 70 anticommutes 
with _ff/f[74, 0] and acts like an identity matrix in Hq. 
The zero energy states in the presence of the pseudo 
field, therefore, differ from the one in its absenc e by 
a factor exp[p{x)], where b = eijdittj = d'^p{x) [37| . 
However, two zero modes are identical for each spin 
component and have equal characteristic lengths. The 
explicit form of the zero modes without the pseudo fields 
can be computed simply by setting 6 = in the previous 
exercises. 

The exclusive members of the CZ(3) algebra of the 
masses are N = {a^jo), three components of the 
Neel vector, whereas, that of the U(l) part is /2 €5 70 : 
charge-density-wave order. The Zeeman coupling 
A ((73 €5 ^4) is proportional to N3 within "Ho, thus may 
develop its finite expectation value by splitting Hq. On 
the other hand, a charge-density wave order can be 
developed inside the core of the vortex if it is charged. 
Contrary to the previous examples, all four orders 
anti-commute with the Hamiltonian Hk[A, a] even when 
an axial field is present. One can choose T = I2® loi 
for all four order parameters. Hence, their expectation 
values can be computed only from the zero energy states. 



VIII. SUMMARY AND DISCUSSIONS 

To summarize we find that in the presence of real and 
pseudo magnetic fields, all the U(l) symmetric orders 
in graphene can host two zero modes with underlying 
vortex defects. However, two zero modes can have 
either one or two characteristic lengths, depending on 
the nature of the underlying insulating orders. We 
show that, if the underlying order parameter couples 
two inequivalent Dirac points, two zero modes have the 
same characteristic length, whereas, with translational 
symmetric order parameters, zero modes enjoy different 
length scales. The existence of two zero modes leads 
to additional competing order parameters in the core 
of the vortex. These order parameters together close 



a CZ(3) X U{\) algebra, and their expectation values 
always arise only from zero energy sub-space. When the 
mid-gap states have different length scales, e.g., with 
underlying AF and QSH orders, charge is continuous 
and oscillatory] the overall neutrality of the system 
is preserved only far from the vortex core. Mid-gap 
states with identical length scales do not lead to any 
excess charge anywhere in the space. The excess 
local charge can be measured by a scanning tunnel 
microscope (STM) probe. It can, therefore, also serve 
the purpose of litmus test to determine the possible 
broken translational symmetry in the system. Zero 
modes with different scales can be delocalized from each 
other by tuning the relative strength of two magnetic 
fields, thereby keeping a single state in the core. For 
example, when B ^ b, ^'i^o is highly localized near 
the vortex core. However, ^'2,0 ~ 6xp (— |A(i)|di) 
experiences effectively zero field, therefore, exists even 
beyond the vortex core. Moreover, ii B,b ^ Ag/i?, 
q+{r) is localized very close to the vortex core, and 
p+(r) {B + b)R/Aoq+{r) > [111- A finite 

magnetization, therefore, persists beyond the vortex 
core and points opposite to that in the core. Net 
magnetization is still zero. 

Besides the insulating orders, one can also consider 
the underlying vortex of various superconducting orders 
(SCOs) in graphene. Aside from the Kekule super- 
conductors, the rests, e.g., s-wave or f-wave SCOs, 
couple two inequivalent Dirac points and may not host 
Majorana modes [1, [3l|. However, a Kekule SCO with 
topological defects, e.g., vortex or half- vortex can lead to 
Majorana modes Those Majorana modes are valley 
polarized. Hence, in the presence of real and pseudo 
magnetic fields, the Majorana modes experience different 
effective magnetic fields. Therefore, by adjusting the 
relative strength of the fields, one can localize one 
of the Majorana modes within the vortex core, while 
the other one can be pushed outside. Binding of a 
single Majorana mode can also be useful for quantum 
computations. Otherwise, the existence of zero modes 
in the presence of vortices can lead to KT scaling of 
the longitudinal resistivity (Rxx) in a neutral graphene 
(ly ~ 0), subject to real and pseudo magnetic fields 
[15| . Zero modes with different length scales can also be 
found in the spectrum of a birefringent Dirac fermion 
with a kinked staggered chemical potential [s^. In a 
recent paper [i^, authors considered dual vortices: a 
vortex in the birefringent parameter as well as the stan- 
dard mass vortex. Under that circumstance, vortex and 
anti- vortex zero modes decay with different length scales. 
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